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Abstract. Hyperon resonances are becoming an extremely useful tool allowing the study of the 
properties of hadronic fireballs made in heavy ion collisions. Their yield, compared to stable 
particles with the same quark composition, depends on hadronization conditions. The resonance's 
short lifetime makes them ideal probes of the fireball chemical freeze-out mechanisms. An analysis 
of resonance abundance in heavy ion collisions should be capable of distinguishing between possible 
hadronization scenarios, in particular between sudden and gradual hadronization. In this paper, 
^vq ■ we review the existing SPS and RHIC experimental data on resonance production in heavy ion 

' collisions, and discuss in terms of both thermal and microscopic models the yields of the two 

observed resonances, K* and A( 1520). We show how freeze-out properties, namely chemical freeze- 

■ out temperature and the lifetime of the interacting hadron phase which follows, can be related to 
resonance yields. Finally, we apply these methods to SPS and RHIC measurements, discuss the 
significance and interpretations of our findings, and suggest further measurements which may help 

■ in clarifying existing ambiguities. 

o 

*cL: INTRODUCTION AND MOTIVATION 

PS. 

The experimental measurement of short-lived hadron resonances can potentially be 
very useful in clarifying some of the least understood aspects of heavy ion collisions. 
In general, evolution of a hot hadronic system proceeds according to Fig. 1. When 
mesons and baryons emerge from a pre-hadronic state, presumably quark gluon plasma, 
their abundances are expected to be fixed by hadronization temperature and chemical 
fugacities. This stage of fireball evolution is commonly known as chemical freeze-out. 
After initial hadronization, the system may evolve as an interacting hadron gas. At a 
certain point (which can vary according to particle species), thermal freeze-out, where 
hadrons stop interacting, is reached. 

A quantitative understanding of the above picture is crucial for any meaningful anal- 
ysis of the final state particles. Many probes of deconfinement are most sensitive when 
the dense hadron matter fireball breakup is sudden and re-interaction time short or non- 
existent. Final state particles could, however also emerge remembering relatively little 
about their primordial source, having been subject to re-scattering in purely hadronic 
gas phase. In fact, theoretical arguments have been advanced in support of both a sud- 
den reaction picture [1,2] and a long re-interaction timescale [3]. Both pictures have 
been applied to phenomenological fits of hyperon abundances and distributions [4, 5]. 

It is apparent that hyperon resonances can be crucial in resolving this ambiguity: 
their initial abundance, compared to stable particles with the same quark composition, 
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FIGURE 1. Stages of the space-time evolution of a heavy ion collision. At a certain moment in proper 
time (known as chemical freeze-out), hadrons emerge. The system then evolves as an interacting hadron 
gas, until thermal freeze-out, the point at which all elastic interactions cease as well. 



depends primarily to the temperature at hadron formation. However, the observed abun- 
dance will be potentially quite different, and will strongly depend on the lifetime of the 
interacting hadron gas phase. Resonances can only be observed via invariant mass recon- 
struction, and their short lifetime means that they can decay within the interacting hadron 
gas (Fig. 2). In this case, the decay products can undergo re-scattering, and emerge from 
the fireball with no memory about the parent resonance. Thus, the observable resonance 
abundance is sensitive to precisely those parameters needed to distinguish between the 
sudden and gradual freeze-out models. 

This paper offers a pedagogical description of what we know about resonance produc- 
tion in heavy ion collisions. We start with a review of presently available experimental 
data, and the open questions which arise. We then proceed to describe how to calculate 
the initial resonance abundance and the effect of re-scattering in terms of the hadroniza- 
tion temperature and the lifetime of the interacting hadron gas. We show how these two 
parameters can be extracted from the experimental observations. Finally, we discuss 
possible answers to experimental challenges raised in the first section, and suggest ways 
by which these questions could be resolved by further measurements. This survey in- 
corporates recently published experimental measurements [6, 7, 8, 9, 10] and theoretical 
results [13, 14, 15]. 



Strange Hadron Resonances: Freeze-Out Probes in Heavy-Ion Collisions February 7, 2008 2 





m=ZE 2 -(Ep) 2 ' m=EE 2 -(Ip) 2 

FIGURE 2. Re-scattering can inhibit resonance reconstruction. Left: in case the resonance decay prod- 
ucts reach the detector without further interactions, their invariant mass distribution should yield a clear 
peak at the resonance mass, with the yield corresponding to resonance abundance. Right: should these 
decay products undergo re-scattering before reaching the detector, the signal may be indistinguishable 
from the background caused by unrelated particle pairs. Generally a weakening by the medium of the 
resonance invariant mass signal must be allowed for. 



EXPERIMENTAL MEASUREMENTS 
Resonance Reconstruction 

Resonances are detected through invariant mass reconstruction of their decay prod- 
ucts. The channels relevant for the here considered K*°, A( 1520) resonances are: 
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FIGURE 3. A(1520) mass plot after mixed event background subtraction with a Breit-Wigner fit from 
NA49 at y/s = 17 GeV, insert plot: Mass plot before mixed event background subtraction [8]. 



Z*°(892) ^k + +K-, 
£*°(892)^7T +K+, 
A(1520) -+p + K~, 



and the proposed E* measurement would be performed via the observation of the decay: 

E*(1385) ^A+n. 

At present, NA49 reported a K* — > K K + signal [6], while STAR published an average 
of if* and W yields [7]. 

Both the STAR [16] and the NA49 [17] experiments perform this reconstruction based 
on data collected in their large volume Time Projection Chambers (TPC) detectors. 
Charged particle momenta are obtained by a measurement of their trajectories in a 
uniform magnetic field, while particle identification is done evaluating the rate of energy 
loss dE/dx in the TPC [18]. A centrality trigger selects the 14% (STAR) and 5% 
(NA49) most central inelastic interactions. The decay daughter candidates are selected 
via their momenta and dE/dx. The resonance signal is obtained by the invariant mass 
reconstruction of each pair combination and the subtraction of a mixed event background 
estimated by combining candidates from different events ( Figure 3). The multiplicity 
is obtained after applying correction factors for acceptance, particle identification cut, 
reconstruction efficiency and branching ratio. 
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TABLE 1. Resonance mass and width 



collision 


Energy 


particle 


mass 


Cmass 


width 


"'width 


reference 




GeV 




MeV/c 2 


MeV/c 2 


MeV/c 2 


MeV/c 2 




Pb+Pb 


17.2 


A(1520) 


1518.1 


2.0 


22.7 


6.5 


[8] 


Au+Au 


130 


K*° 


893.0 


3.0 


58.0 


15.0 


[10] 


Au+Au 


130 




896.0 


4.0 


63.0 


11.0 


[10] 


PDG 




A(1520) 


1519.5 


1.0 


15.6 


0.6 


[11] 


PDG 




K*° and F° 


896.1 


0.27 


50.7 


15.0 


[11] 



Determination of Width and Mass 

The measured width and mass of A(1520) and K* in heavy ion collisions is in 
agreement within the statistical errors with the values from the Particle Data Group 
(table 1). 

Within the statistical error there are no width broadening or mass shifts observed. 
While there is no calculation available that estimates the width and mass profile of a 
resonance by taking the medium density evolution during the expansion into account, 
shifts in width and mass would be expected in a lengthy hadron gas phase [19]. A 
prediction based on relativistic chiral SU(3) dynamic calculations [20, 21] gives a 100 
MeV broadening of the A(1520) resonance and a mass shift of about 100 MeV to lower 
masses in medium at a density of p = 0.17 fm~ 3 . 

The fact that the A(1520) NA49 signal from Pb+Pb reactions shows no width broad- 
ening within the errors, indicates that only the decay products coming from the vacuum 
A(1520) decay are observed. Any medium modified A(1520) resonances must thus de- 
cay rapidly within the medium, with decay products rescattering, or, alternatively no 
such A(1520)-medium interaction is present. To distinguish these alternatives one needs 
to consider the yields of A(1520). 



Particle Multiplicities 

Experimental results for K* and A(1520) resonance production in heavy ion colli- 
sions are available from NA49 experiment in central Pb+Pb and p+p at s/s =17 GeV 
and from STAR experiment in central Au+Au y/s^ =130 GeV collision energy and 
are shown in table 2. 

Comparing A(1520) — > pK yields at different experiments, we see a marked decrease 
from elementary p+p to Pb+Pb collisions at 17.2 GeV, which is even more pronounced 
in the A(1520)/A ratio (figure 4 left and right). This suggests that heavy ion collisions 
are not merely superpositions of elementary p+p collisions. Here we can ask the first 
question. "Why is the A(1520) signal from the invariant mass reconstruction in the K-p 
decay channel actually suppressed in heavy ion collisions, whereas all other hyperons 
are enhanced?" Moreover, the K* over % production at the same energy indicates no such 
signal decrease from p+p to Pb+Pb collisions, shown in figure 5. The second question is 
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TABLE 2. Resonance yields 



Energy 



collision 


GeV 


particle 


yield 


reference 


p+p 


17.2 


A(1520) 


0.0121 ± 0.003 


[6,8] 


p+Pb 


17.2 


A(1520) 


0.0072 ± 0.003 


[6] 


Pb+Pb 


17.2 


A(1520) 


1.45 ±0.40 


[6,8] 


Pb+Pb 


17.2 


K* 5-10% most central 


12.88 ±4.53 


[6] 


Au+Au 


130 


K* (dN/dy) y<|0 5| 


10.0 ±0.8 


[10,7] 


Au+Au 


130 


A(1520) 


<4.2* 


[9] 



* at 95% confidence level 
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FIGURE 4. Left: A(1520)/7T+ as a function of participants for p+p, p+Pb and Pb+Pb collisions at 
158 AGeV [6]. Right: (A(1520))/A ratio as a function of number of participants. Data points are p+p 
and central Pb+Pb collisions from NA49. 



now "Why is no comparable suppression observed in the K*7" A hint for the signal loss 
could be due to secondary interactions is the measurement of the A(1520) production 
in p+Pb collisions which seems to be decreasing. However the errors are to large to say 
this conclusively (figure 4 left). 

For the collision energy of y/s^ =130 GeV there is only an upper limit for the 
A(1520), namely 4.2 at 95% confidence level (2a). The expected multiplicity from 
extrapolated elementary p+p reactions including an addition factor of 2 for strangeness 
enhancement (taken from A at SPS [22]) is ~ 7.7. The upper limit estimate indicates 
that we see at RHIC energies the same trend of signal loss as at SPS energies. The K*/K 
ratio at RHIC energies fits in within errors to the ratios from elementary p+p and e + +e 
collisions figure 6, that leads us to the conclusion that there is no indication of K* signal 
loss related to K yield. Our conclusions from the results from SPS and RHIC energies 
are in agreement. 

A recent UrQMD calculation [33] predicts for the K* a signal loss of 66 % and for the 
A(1520) 50 % at SPS energies due to re-scattering of the decay products (figure 7). At 
RHIC energies the signal loss is on the order of 55% for K* and 30 % for A(1520) [34]. 
The K* yield at A+A collisions for SPS and RHIC energies shows little or no indication 
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FIGURE 5. K*(892)/7T as a function of participants for p+p and Pb+Pb collisions at 158 AGeV. 
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FIGURE 6. K*°/K measured in different colliding systems at different energies. For elementary e + e 
and pp collisions in comparison [23, 24, 25, 26, 27, 28, 30, 31, 32] with the heavy ion collision Au+Au at 
RHIC energies = 130 GeV [7]. 



of signal loss which was expected. This means that the mechanism of re-scattering alone 
can not explain the signal loss of the A(1520) with respect to the K* production. 

As table 3 shows, the agreement of the data with thermal model predictions presents a 
considerable problem as well. At SPS energies, thermal model calculations over-predict 
theA(1520) yield by a factor of two [36], while the upper limit measurement of A( 15 20) 
from STAR is 30% lower than the predicted value. However, published thermal model- 
calculations significantly under-predict the K* (although, as will be seen in the next 
section, the observed K* is consistant with thermal production at a lower freezeout 
temperature). 
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TABLE 3. Particle yields from thermal model predictions 





Energy 




predicted 




collision 


GeV 


particle 


yield 


reference 


Pb+Pb 


17.2 


A(1520) 


3.48 


[35] 


Pb+Pb 


17.2 


A(1520) 


5.20 


[3] 


Au+Au 


130 


K*/hT 


0.037 


[12] 




MODELING RESONANCES IN HEAVY ION COLLISIONS 
Direct production at hadronization 

We assume that at hadronization, a volume element will be at thermal and chemical 
equilibrium in its local rest frame. This means hadrons produced directly from a medium 
at temperature T (much lower than the particle mass, as in all cases considered here) fill 
the available phase space according to the relativistic Boltzmann distribution: 

dNocgd 3 pf(X,E,T), (1) 
f{X,EJ)=Xe- E l T . (2) 

Here g is the statistical degeneracy and X = e^/ T is the hadron's fugacity. Since a 
hadron is a composite state of flavored quarks, it is natural to parametrize its fugacity 
as a product of the flavor fugacities of each constituent. As a further step, each quark's 
fugacity can be decomposed into a term X- controlling the net flavor (q — q) and an 
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occupancy number y. which determines the number of quark- antiquark pairs. Therefore, 









q=l 




= v. 




= y T 



(3) 

(4) 
(5) 

and the precise values of X- v y i can be determined from global particle abundance fits [2]. 

We can now change to rapidity (y) and transverse mass (m r ) coordinate system, 
where transverse and longitudinal directions are defined with respect to the beam. (The 
advantage of this system is that m T is invariant between the lab and the center of mass 
frames, while the rapidity transforms additively). 



dN = gm T cosh(y)dm T dyf(X,m T cosh(y),T) (6) 



y = hn(^) (8) 



2 \E-p L 

To get particle ratios in a certain region of phase space collective expansion (flow) as well 
as freeze-out geometry has to be included. The most general way to do this is to introduce 
a space-time hypersurface (a 3-D surface in 4-D space-time) dZ^, which will transform 
as a 4-vector. This vector can either be time-like, corresponding to a freeze-out surface 
evolving in time, or space-like, corresponding to a freeze-out occurring simultaneously 
across the fireball volume. The total number of emitted particles will then be 

N l = JdZ^x)f(x). (9) 

We can describe the current 4-vector in terms of the hadron density in the volume 
element at rest, (see Eq. 1) 

/ = n rest ^ = /(A,£ local ,T)^, (10) 

where * is the energy in the local rest frame with respect to flow, and is the 
volume element's 4-velocity. 

In the limit of a homogeneous fireball (Temperature and chemical potential constant 
in space) and full phase space, the total number of produced particles is, as expected, 
independent of flow and hadronization geometry, since in this case Eq. (9) becomes 

N = j dL ll (x)f{x) = J d^(x)u^x) J d 3 pf(?i,E local ,T) = vJ d 3 pf(X,E local ,T). 

(ID 

We can now find N analytically 

N = VkJ d 3 pe-( p2+m2 V T = ^Xm 2 TK 2 {m/T) (12) 
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remembering the properties of Bessel functions 



Ux) = p-x~ n / dz(z 2 + x 2 )-h 2n e-^ 
{In ) ! Jo 



(13) 



If the fireball is not homogeneous, or we want to only calculate the yield within the 
experimental acceptance limits, we need to use a more general result. Putting Eqs. (9) 
and (10) together one gets the relativistically invariant Cooper- Frye formula [40] 

= _ 1 f N = [ dL^f(X,p^,T)G(L^), (14) 
aPp 27tm T aym T J 

where the step function 0(JL^p^) has been introduced to eliminate unphysical emission 
in the direction opposite to expansion [41, 42]. Eq. (14) can be integrated to yield the 
total particle number over an arbitrary phase space region 

Many of the parameters used in these formulae can be eliminated by judicious choice 
of observables. For instance, the ratio of a hyperon resonance to its ground state is 
independent of the chemical potential, since all chemical potential terms cancel out. 

Suitable ratios include E*/A, A(1520)/A, K*°/K+, K^/K'. A particular case is the 

K°/K° system: While the decays long before it has a chance of oscillate, the 

experimentally observed K s is a superposition of K° and K°. Suitable average are 

K^/K- (measured by NA49 [6] ) K*°/K+ and (K* +W)/(K+ +K~) (measured by 
STAR [7]). 

Moreover, we found that for these ratios (more generally for ratios of particles with 
comparable masses), effects due to flow and surface geometry cancel out to a very good 
approximation and Eq. (12) gives the produced particle ratios to a good accuracy. 



Feed down from resonance decays 

Many of the hyperons considered here can be produced either directly at hadronization 
or in the decay of a resonance produced earlier. For instance, 

A(1520) A(1520) 

A total A+(E°^Ay) + (E* ^Att) + (E* ^YPk^Akj)' 

Unless a resonance will be directly reconstructed by the experiment, its decay products 
will be a component in the hyperon abundances (albeit with a subtly different momentum 
distribution). This contribution, therefore, will need to be taken into account through 
calculation. As in the case of direct production, all kinematic effects are integrated out 
in a calculation of the total yield, and one simply sums each term in the resonance 
contribution as a term of the form given in Eq. (12). If a finite acceptance in phase space 
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needs to be taken into account, however, a distribution of decay products will need to be 
computed from the statistical distribution of resonances. 
We assume that, in a decay of the form, 

*->l+2+..., (17) 

any dynamical effects in the decay (the S-matrix) average out over a statistical sample 
of many resonances. In other words, in the rest frame comoving with the "average" 
resonance, the distribution of the decay products will be isotropic. The rate of particles 
of type 1 (as in Eq. (17)) produced with momentum p* in the frame at rest w.r.t. the 
resonance will then simply be given by the Lorenz-Invariant phase space factor of a 
particle of mass M x and momentum p* within a system with center of mass energy 
equal to the resonance mass M R . 

d 3 N { " " ^ 3 "* N 



d 3 p] 



b I fl^S S (LPi -P^CLEt-Et ~M R ), (18) 



where b is the branching ratio of the considered decay channel. All that is left is to 
change coordinates from the resonance's rest frame (p*, E*) to the lab frame (p,E). 

If more than two bodies are considered, this calculation becomes more involved [45]. 
For the general N-body case, it is better left to Monte Carlo methods [46]. In the case of 
the 2-body decay, the situation is greatly simplified by the fact that p*, E* are fixed by 
energy-momentum conservation to, 

Et = ^ f (M 2 R -m 2 l -m 2 2 ), (19) 

R 



pX = -pt = ^Ef-m\. (20) 
Putting the constraints in Eq. (19) into Eq. (18) one gets, after some algebra [44] 

dN b f Y + f M r+ , d 2 N R 

L ^-ttt^t, (21) 



dm 2 ,dy 1 Anp\ Jy- Jm t - dM 2 , R dY R 

J = M * (22) 



P TRPT\ ~ ( M R E R - M TR m T\ COSh AY) 2 

A7 = Y R -y v (23) 

J is simply the Jacobian of the transformation from the resonance rest frame to the lab 
frame, and the limits of the kinematically allowed integration region are: 



l / Pi 



! 



Y ± =y x ±sinh 



E R m T1 cosh(Ay) ± p T1 J p* 2 - mj, l sinh 2 (Ay) 

= M R =-i . 

Ml m — 



m; 



, 1 sinh 2 (Ay)+m^ 
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TABLE 4. Resonances contributing to A and K production, 
with their degeneracies, rest-frame momentum (p*) and possi- 
bility for experimental reconstruction 



g 


Reaction 


F 
MeV 


branching 


visible? 


«4 


E*°(1385)^lV 


127 


«4% 


No 


8 


E* ± (1385) — * Arc* 


208 


88% 


Yes 


4 


£*°(1385)^A7r 


208 


88% 


No 


2 


E° — > Ay 


74 


100% 


No 


4 


A(1520)^NF 


244 


45% 


Yes 




3 


AT* (892) -fAT+Tr 


291 


67% 


Yes 



We now apply the methods outlined here to calculate ratios at hadronization for a 
range of experimentally observable resonances, both at full and central rapidity (|y| < 
0.5). The total ratio was calculated using the analytical formula in Eq. (12), while for 
mid-rapidity equations in Eq. (15) and (21) were used. In the midrapidity calculation we 
used a uniform spherically symmetric expanding fireball, with a hadronization surface 
used earlier by fits in [5] and [44], 

dL fl = (l-'^e r )d s ? t (24) 
r or 

where e r is the unit vector in the radial direction and dt J dr is a constant. We mentioned 
earlier that flow and hadronization surface details cancel out to a very good approxima- 
tion when hadron masses are comparable. In fact, when we varied dt^/ dr, the results did 
not change by more than a percentage point. Table 4 summarizes the decay processes 
considered in our analysis and their parameters (Clebsh- Gordon coefficients have been 
used to estimate decays such as (N*° ^N + n~)/ (N*° — > A^ )). 

In Fig. 8 we show the relative thermal production ratios at chemical freeze-out over 
the entire spectrum of rapidity and m T (solid lines) and central rapidity range (dashed 
lines). The sensitivity of resonance yields to hadronization temperature is apparent for 
all resonances under consideration. In particular, the E* emerges as a very promising 
candidate for further study. For example, at the lowest current estimates (T ~ 100 MeV) 
of the final break up temperature in 158A GeV SPS collisions 33% of A are arising 
from primary £*, the percentage rises to slightly more than 50% if chemical freeze-out 
occurs at T = 190 MeV. Both SPS and STAR detectors are capable to measure hyperon 
and resonance yields well within the precision required to distinguish between the two 
limiting cases. 

as we discussed earlier, the experimentally observed K* yield is compatible with the 
thermally produced ratio for a range of hadronization temperatures, but the A(1520) 
seems very suppressed. To discuss this further, an estimation of the effect of rescattering 
on resonance abundances is necessary. 
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FIGURE 8. Temperature dependence of ratios of £*, K* and A( 1520) to the total number of observed 
K A s and E s. Branching ratios are included. Dashed lines show the result for a measurement at central 
rapidity Ay = ±0.5. The experimental measurements (horizontal lines, not E* /E) included in the diagram 
were presented in the first part of this paper. For the E yields in the diagram, chemical potentials were 
taken from [2]. It should be noted that the /if* //if ratio is actually K* /K~ for the NA49 measurement, and 
an average (K* +K*)/ (K+ + K ) for STAR. 



RESCATTERING 

As explained in the introduction, direct observation of resonances relies on invariant 
mass reconstruction. Therefore, to calculate the observed resonance abundances, rescat- 
tering after hadronization will have to be taken into account. This can be looked at within 
a microscopic model of hadronic matter such as UrQMD. As we have shown earlier, 
UrQMD alone did not explain consistently the experimental data if a long re-interaction 
period and initial thermal model yields were assumed. However, as yet there was no de- 
tailed study of the dependence of the resonance yields on the lifetime of the interacting 
hadron gas phase (the interval between chemical and thermal freeze-out) as well as the 
chemical freeze-out temperature. 

Here, we present such a study using a "back of the envelope" model which neverthe- 
less seems to provide an acceptable qualitative description of the propagation of reso- 
nances and their decay products through opaque matter. 

We first note in Fig. 8, that the relative E*/S signal is remarkably independent (within 
5%) of Temperature. This is because the S*(1530) contribution cancels nearly exactly 
the thermal suppression of the S originating in the S — E* mass difference. This effect 
could be used for a direct estimate of the E* lost through rescattering, even without 
knowledge of the freeze-out temperature, should the chemical parameters (X l and j\) 
be independently known. A simple test of sudden hadronization model consists in 
measurement of the ratio E*/S. If it is significantly smaller than unity, we should expect 
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a re-equilibration mechanism to be present. Otherwise sudden hadronization probably 
applies, since E* emerge from chemical freeze-out without undergoing interactions. 

We can, however, go further and use the suppression of the considered resonances as 
a tool capable of yielding a quantitative estimate of the lifetime of the interacting hadron 
gas phase. We consider the decay of a generic resonance Y* 



Yk. 



(25) 



in a gas of pions and nucleons. We shall assume that one interaction of either Y or % 
is sufficient for that resonance to be undetectable, and that the decay products travel 
through the medium with speed \ i (where i can mean either Y or %). The interaction 
probability is proportional to \- v the interaction cross-section of the decay product with 
each particle in the hadronic medium (cr^vj, where j can refer to either pions, Kaons, 
nucleons or antinucleons. Note that the cross-section itself depends, in a generally 
complicated way, on the incident momentum, and hence on v 4 ), and the particle density 

in the fireball pj. pj is increased by a factor y { = 1/ y^l — v? due to Lorentz-contraction, 

and decreases as time passes because of the fireball's collective expansion (parametrized 
by the flow velocity v flow , assumed to be of the order of the relativistic sound speed 

cj \/3.) The time dependence of the densities will therefore be, 



R 



R + \ 



flow' 



(26) 



and p ., the density of j at hadronization, can be calculated from the chemical freeze- 
out temperatures and chemical potentials. Putting everything together, the rescattering 
reaction rate is 



P i = E a ^( V i)P ^ + ^dPoK + <%( V i)P(W + %( V i)P, 



0N J 



If we use the average, 

L cT ( v i) v iH 



.< a >< vv. >=< a > — , 



(rv)i 



R 



R + \ 



flow' 



(27) 



(28) 



(where p and m are the resonance's momentum and mass) Eq. (27) becomes, 



Pok + (^k) Pok + <<%> Pon + (%) P 0N 



1- 









(29) 



Neglecting in-medium resonance regeneration and particle escape from the fireball, the 
population equation describing the scattering loss abundance (N { ) is: 



dt 



1,2 



dN, 



dt 



— N, 



(30) 
(31) 
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TABLE 5. Scattering model parameters 



'nN 



(mb) 



24 



20 40 



20 



24 



NN 

50 



35 MeV 



x A(1520) 

15.6 MeV 



£*°(892) 

50 MeV 



escape rate (fm 3 ) 



negligible 



0.5 



R(fm) 



8-145/r [MeV] 



220 MeV 



The required nucleon and antinucleon density at hadronization, p QN is obtained through 
Eq. (12) 

Po » = (2^ W(A ^ )37 ^ ( ? ) ' (32) 

We consider the nucleons to have a mass of ~1 GeV, and a degeneracy of six, to take the 
p,n and the thermally suppressed but higher degeneracy A contributions into account, 
the pion density at hadronization is computed in the massless particle limit, leading to 

Po, = ^ 3 , (33) 

The model presented here is remarkably insensitive to the individual cross-sectional 
areas. The values we used in the calculation are given in table 5, but order-of-magnitude 
variations of the more uncertain cross-sections did not produce variations of more than 
30%. Similarly, the value of the initial fireball radius R Q (which is constrained by the 
entropy per baryon) does not significantly affect the final ratios. This reassures us that 
had we used a more exact approach than the approximations in Eq. (28), the qualitative 
features of our model would not have changed. The results, however, exhibit a very 
strong dependence on both the Temperature (which fixes the initial resonance yield as 
well as the hadron density of the fireball) and fireball lifetime (in a short-lived fireball 
not many resonances decay, so their products do not get a chance to rescatter). Figure 9 
shows the dependence of the A(1520) /A, E*/A and K*°(S92) /K on the temperature and 
lifetime of the interacting phase. It is clear that, given a determination of the respective 
signals to a reasonable precision, a qualitative distinction between the high temperature 
chemical freeze-out scenario followed by a rescattering phase and the low temperature 
sudden hadronization scenario can be made. We also note that despite the shorter lifetime 
of the E* and higher pion interaction cross section, more E* decay products should 
be reconstructible than in the A(1520) case, at all but the highest temperatures under 
consideration. This reinforces our proposal that the E* is a very good candidate for 
further measurement. 

Diagrams such as those in Fig. 9 still contain an ambiguity between temperature and 
lifetime of the interacting hadron gas phase. A low observed ratio can either mean a low 
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T [MeV] T [MeV] 

FIGURE 9. Produced (dashed line) and observable (solid lines) ratios A(1520)/(total A),£7(total 
A), K /K and S/(total A). The solid lines correspond to evolution after chemical freeze-out of 
1,2,3,4,5,7,10,15,20 fm/c, respectively. The values at time zero (chemical freezeout) were taken from 
Fig. 8. See Fig. 8 caption for the meaning of K,K*. 



freeze-out temperature or a lot of rescattering in a long re-interaction phase. However, 
this ambiguity can be resolved by looking at a selection of resonances, with different 
masses and lifetimes. Fig. 10 shows how the initial temperature and the lifetime of the re- 
interaction phase decouple when two resonance ratios are measured simultaneously. A 
data point on diagrams such as those in Fig. 10 is enough to measure both the hadroniza- 
tion temperature and to distinguish between the sudden freeze-out scenario and a long 
re-interaction phase. The plots in Fig. 10 can also be used as consistency checks for 
the model: For example, the near independence of E*/S on temperature means that the 
equal temperature lines in the E*/S vs E*/A diagram are nearly insensitive to the details 
of the rescattering model. Moreover the mass differences and lifetimes of the E*, K* 
combine in such a way as to make the E*/(all A) vs K*° (S92) / (all K-) diagram fold into 
a very narrow band. Any serious shortcoming within our rescattering model would be 
revealed if the observed particle ratios stray from this band. 
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FIGURE 10. How temperature and fireball lifetime decouple when two resonances of different masses 
and widths are compared. A point on any of the above diagrams is potentially sufficient to fix both of 
these quantities. The experimental measurements are discussed in the first section. See Fig. 8 caption for 
the meaning of K,K* . 



DISCUSSION 

The first experimental results on short-lived resonances have raised more questions 
than answers. The low A(1520) multiplicity measured by NA49, together with its non- 
detection at STAR, can not be understood with a thermal model exclusively. A A(1520) 
suppression factor of at least 2 is needed to reproduce the measured multiplicity. How- 
ever, this suppression does not appear to result from simple scattering with in-medium 
particles, described by models such as UrQMD. Such suppression should manifest it- 
self much more strongly in case of the K* , which has a much shorter lifetime. On the 
contrary, no comparable depletion has been found in the case of the K* either at SPS or 
RHIC energies. Moreover, no broadening of the resonance width, which should accom- 
pany a strongly interacting medium, has been observed. If these results are incorporated 
in the diagrams of Fig. 10 , the apparent result is sudden freeze-out with an unphys- 
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FIGURE 11. K*/K vs A(1520)/A assuming half of the A(1520) s are suppressed in-medium. See Fig. 
8 caption for the meaning of K,K*. 



ically low (~ 100 MeV) temperature. In analysing these results it should be kept in 
mind that the A(1520) is a very unusual particle: Unlike most other hadronic resonances 
(E*, K* , A, p ecc), its high spin is due not to valence quark spin configuration but to the 
fact that the A(1520) s valence quarks are believed to be in an L=l state. The greater 
space separation of L=l wave-functions means that the A(1520) is especially suscepti- 
ble to in-medium modifications. As Fig. 11 shows, a 50% suppression of the A(1520) 
signal at hadronization would mean the data is perfectly compatible with the sudden 
freeze-out model described in [1] 

Some proposed mechanisms which could lead to a A(1520) suppression (Some sug- 
gested by Berndt Muller at the meeting [49]) are: 

• If hadrons form in a quark coalescence process, the A(1520) would be suppressed 
for the same reason that L=l (p-shell) electrons are never captured by nuclei 
[50]: The coalescence probability is dependent on the wave-function overlap at the 
formed particle's location ~ |\P(;c = 0)| 2 . The overlap of two particles goes down 
for higher L states, and is for odd L=l,3,... states. However, this effect should also 
manifest itself in proton-proton collisions, which serve as basis of the AA-reaction 
yield. We can thus exclude this effect. 

• Similarly, the larger spatial extent of the A(1520) means it is more susceptible to 
in-medium screening of the strong potential. This mechanism invoked for Charmo- 
nium suppression can also suppress states such as the A(1520). 

• In-medium effects can mix the A( 1520) \S = 3/2,1 = 1 /2) with the E* 
\S = 3/2,7 = 3/2) state. Reactions such as 



?r + A(1520) — > E* 



(34) 
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can deplete the A(1520) states without a significant effect on the (much more 
abundant) E* population [14]. This is analogous to the phenomenon of quenching, 
well known in atomic resonances [51]: Resonance radiation from many gases is 
known to be suppressed by in-medium collisions, which prevent photon emission 
by converting the excitation energy of the resonance into kinetic energy. 

• Spin-orbit interactions due to in-medium chromomagnetic fields can also mix the 
A(1520) with the E*. 

• Isospin conservation forbids the decay A(1520) — > An but not A(1520) — ► Ac, 
which is forbidden in vacuum since m A + m a > fn A ^ 520 y At a temperature high 
enough for (partial) chiral symmetry restoration to take place, this may no longer be 
true, and the A( 1 520) — > Aa decay would become possible. In the same way, partial 
restoration of chiral symmetry may be suppressing the K* — > Kn decay channel. 

• Finally, in a long re-interaction phase there is a possibility that the K* will be re- 
generated. Processes such as 

Kk -> K* -> Kk 

will lead to an enhancement of observed K*. Processes which regenerate the 
A(1520), such as 

pi^^A(1520) -> pK 
are considerably less likely, due to the suppression of higher (L > 0) partial waves. 

The last effect can be properly taken into account by a more elaborate rescattering 
model, which takes all of the microscopic reaction dynamics into account. The other 
mechanisms are rather complex and uncertain, though for the RHIC system, potential 
screening, spin-orbit mixing and o mass change can be explored in a \l B = lattice 
calculation. However, these effects also depend on the unusual characteristics of the 
A(1520) resonance. It seems to us that all these consideration suggest that other reso- 
nances should be studied to constrain the fireball freeze-out properties. 

The non-suppression of the K* makes it likely that other hadron resonances can also 
be detectable. The diagrams shown in the previous section imply that the E* should be 
abundantly produced, and it has many characteristics which would make it a logical next 
step in the study of resonances produced in heavy ion collisions. Other proposed targets 
of observation are the p and even the E° — > Ay. 

The sudden freeze-out scenario can be thoroughly tested by the detection of non- 
strange baryon resonances (N* ( 1440) and A( 1230)). The larger widths of these particles, 
as well as the larger scattering cross-section in hadronic matter of their decay products, 
mean that the very detectability of these particles in terms of an invariant mass analysis 
would be strong evidence of a very fast hadronization process. 
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